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Abstract. The behaviour of solitons in integrable theories is strongly constrained by 
the integrability of the theory; i.e. by the existence of an infinite number of conserved 
quantities which these theories are known to possess. One usually expects the scattering 
of solitons in such theories to be rather simple, i.e. trivial. By contrast, in this paper 
we generate new soliton solutions for the planar integrable chiral model whose scattering 
properties are highly nontrivial; more precisely, in head-on collisions of N indistinguish- 
able solitons the scattering angle (of the emerging structures relative to the incoming 
ones) is 7r/N. 

We also generate soliton-antisoliton solutions with elastic scattering; in particular, a head- 
on collision of a soliton and an antisoliton resulting in 90° scattering. 



I. Introduction. 



This paper studies certain exact soliton solutions of an integrable system. Before any 
detailed discussion, and to avoid confusion later on, it is worthwhile clearing up a small 
point of terminology: the word solitons was introduced by mathematicians to describe 
lumps of energy which were stable to perturbations and did not change either velocity or 
shape when colliding with each other. However, in recent literature all sorts of localized 
energy configurations have been called solitons. We shall go along this looser definition. 
By a soliton we shall mean a lump of energy that moves but we shall not imply stability 
of the shape or the velocity or a simple behaviour in collision. 

An interesting problem is to look at the scattering properties of two or more solitons 
colliding. In some known systems with nontrivial topology, the collision of two solitons 
is inelastic (some radiation is emitted) and nontrivial (a head-on collision results in 90° 
scattering); all this has been observed analytically |], @] and numerically @-|§. One 
can construct explicit time-dependent solutions only in very special, so-called integrable 
models. Usually in these models extended objects interact trivially, in the sense that they 
pass through each other with no lasting change in velocity or shape (i.e., they behave as 
genuine solitons). Some examples in (2+1) dimensions are the Kadomtsev-Petviashvili 
equation J7J and the integrable chiral model 0. The last system is the subject of this 
paper and will be described below. Until now, nontrivial scattering of solitons occurs 
mostly in nonintegrable systems which is far from simple. The question that arises is 
whether this type of scattering can occur in integrable models too. There are some 
limited examples of integrable systems where soliton dynamics can be nontrivial. In (1+1) 
dimensions there many models which possess nontrivial soliton-like solutions (cf. ||); like 



the boomeron solutions [iQl, which are solitons with time dependent velocities. In (2+1) 



dimensions there are the dromion solutions [11] of the Davey-Stewartson equation, which 



decay exponentially in both spatial coordinates and interact in a nontrivial manner |12|] ; 
and the soliton solutions Jl3| of the Kadomtsev-Petviashvili equations, whose scattering 
properties are highly nontrivial. 

In the present work we are going to construct families of soliton solutions for the integrable 
(2+l)-dimensional chiral model and observe the occurrence of different types of behaviour. 
This happens since the solitons in this system have internal degrees of freedom which 



determine their orientation in space; do not affect the initial energy density; and are 
important in understanding the evolution as a whole. Therefore, they can interact either 
trivially or nontrivially, depending on the orientation of these internal parameters and on 
the values of the impact parameter defined as the distance of closest of approach between 
their centres in the absence of interaction. Namely, if two initial soliton-like structures are 
sent towards each other at zero impact parameter, then, as most numerical simulations 
have shown, the outgoing structures emerge at 90°. 

To proceed further let us specify the system. The modified SU(2) chiral model studied 
by Ward || is given by the field equation 



Here J takes values in SU(2) group and is thought of as a 2x2 unitary matrix of functions 
of the spacetime coordinates on R 2+1 : x M = (a; , x 1 , x 2 ) = (t, x, y) with det J = 1. Greek 
letters are spacetime indices, taking values 0, 1, 2, <9 M denotes partial differentiation with 
respect to x^, while J M = d^J . The quantity e a ^ u is the alternating tensor of three indices 
with e 012 = 1. Finally, V a is a unit vector in spacetime. The conformal properties of V a 
determine whether the symmetry group is SO(2) or SO(l,l) (depending on whether V a is 
timelike or spacelike). 

Ward H chooses V a to have the components V a = (0,1,0), the spacelike case, so that 
(0) is a chiral equation with torsion term and has the same conserved energy-momentum 
vector as the chiral field equation. In fact, the corresponding energy density is 



Here tr denotes the matrix trace. It should be emphasized that £ is a positive-defined 
functional of J, and hence a conserved energy exists which is the integral of the energy 
density over the spacelike plane x° = const. The boundary conditions are chosen so that 
the field configuration has finite energy. Consequently, we require that J be everywhere 
smooth and that 



at spatial infinity, with x + iy = re . Here Jo is a constant matrix, and J\ depends only 
on 9 (no time dependence). 




(1) 






(3) 



The ensuing system when V a is i times a timelike vector instead of spacelike has been 
studied in []14 |. Equation (JT]) admits solitons, localized in two dimensions, with trivial 



scattering, i.e. each soliton suffers no change in velocity and no phase shift upon scattering 
0, [14J]. It is the purpose of this paper to construct new soliton solutions for (P, and 



investigate their scattering behaviour. Such solutions are localized along the direction of 
motion; they are not however, of constant size: their height, which corresponds to the 
maximum of the energy density S, is time dependent. 

The rest of the paper is arranged as follows. In the next section we shall briefly discuss 
the integrability properties of (|l|), and write down a family of multisoliton solutions as 
configurations that are the limiting cases of the ones already obtained using the standard 
method of Riemann problem with zeros [|J. In section 3 we construct two families of 
multisoliton solutions with nontrivial scattering; in particular, for the first one we prove 
that in all head-on collisions the N moving structures undergo tt/N scattering. In section 
4 we construct a mixture of soliton-antisoliton solutions, and in section 5 we discuss their 
dynamics and scattering properties. We finish the paper with a short section containing 
our conclusions. 

II. Construction of Soliton Solutions. 

The integrable nature of equation (]l|) means that there is a variety of methods for con- 
structing exact solutions. Together with Riemann problem with zeros ||, both twistor 
techniques |15| and a full inverse scattering formalism [|16j have been applied to the model. 



This section indicates a general method for constructing soliton solutions of the integrable 
chiral model ([!]). The technique is a variation of that in || [17|, following a pioneering 



idea of Zakharov and his collaborators [TLB 



The nonlinear equation ([!]) is integrable in a sense that it may be written as the compat- 
ibility condition for the following linear system 

Lip = (Xd x - d u )tfj = Alp, (A \ 
Mip = (Xd v - d x )tp = Bip. { > 

Here A G C, (w, v, x) are coordinates on R 2+1 with u = (t + y)/2, v = (t — y)/2, A and B 

are 2x2 anti-hermitean trace-free matrices depending only on (u,v,x), and u, v, x) 

is an unimodular 2x2 matrix function satisfying the reality condition 

■0(A, u, v, x) u, v, xy = /, (5) 



where bar denotes complex conjugate, t denotes the complex conjugate transpose matrix 
and I is the 2x2 identity matrix. The system (^) is overdetermined, and in order for a 
solution ip to exist, A and B have to satisfy the integrability conditions, which are 

B X = A V , A x - B u - [A, B] = 0. (6) 

If we put J(u,v,x) = ip{\ = 0, w, x) -1 where ip is a solution of the system (|), we get 
by comparing (§) and (||) that 

A = J' l J v , B = J- l J x . (7) 

Therefore, the integrability condition for (|J) implies that there exists a field J which 
satisfies the equation of motion (P; and moreover, the reality condition on if) ensures that 
J is unitary 

Using the standard method of Riemann problem with zeros, in order to construct multi- 
soliton solution one may assume that the function ip has simple poles in A, or in other 
words must possess the form 



iKa) = j+x; 



k=1 A - /i fc 

where M k are 2x2 matrices independent of the complex parameter A, n is the number of 
solitons, and the complex parameter /i^ determines the velocity of the k-tYi soliton. The 
components of the matrix M k are given in terms of a rational function f k of the complex 
variable u k = x + fi k u + [i k x v. (Roughly speaking, f k {oJh) describes the shape of the fc-th 
soliton). In fact, the matrix M k (cf. H) has the form 

n 

M fe = -E(r1^W, (9) 
i=i 

with T _1 the inverse of 

T kl = jz{nk-^Y l m k a m l a . (10) 

a=l 

Here are holomorphic functions of uj k , given by = (m^mj) = (l,f k ). These 
solitons pass each other without any change of direction or phase shift. Infinite energy 



extended wave solutions |19| may be constructed by taking f k to be an exponential func- 
tion of u k . Such extended wave solutions suffer a phase shift upon scattering, although 
again there is no change in velocity. 



All this assumes that the parameters fik are distinct, and also /ifc 7^ Hi for all k, I. In 
this paper examples are given of two generalizations of these constructions: one involving 
higher-order poles in /ik, and the other where /ifc 7^ fJ>i- 



Let us look at an example in which the function ip has a a double pole in A, and no other 
poles. So we take if) to have the form 

where Rk are 2x2 matrices independent of A. [This hypothesis can be generalized by 
taking the function if) to have a pole of order n in A.] 



It has been proved |T7| that if) given by ([TT|) satisfies the reality condition (|5|) if and only 
if it factorizes as 

™-('-!HliF)('-8^TiF)' (12) 

where q^ are two-dimensional row vectors and ||gfc|| 2 = qk 4 Ok- 

The qt have to satisfy a condition, which amounts to saying the matrices A = (Lip)^' 1 
and B = (Mip)^ -1 are independent of A. One way of obtaining with this property is as 
a limit of the simple-pole case @ with n — 2. The idea is to take a limit — > fi. In order 
to end up with a smooth solution if) for all (u, v, x), it is necessary that f 2(^2) — f 1(^1) 
in this limit. 

In our case, with n = 2, we put //1 — // + e, a*2 — A* — £ an d write f\{oj\) = /(^i), 
f2{^2) = f{ UJ 2), with / being a rational function of one variable. In the limit e — > 0, ip 
has the form (|L2|) , with 

gi = (1 + m 2 )(i, /) + m - i*)(f, -i), , 13) 

<?2 = (!,/)■ 

Here / is a rational function of u; = x + //w + = (it — H~ 2 v) f'(uj), while /'(u;) 

denotes the derivative of f(u>) with respect to its argument. As a result, we have a solution 
J = = O)^ 1 depending on the complex parameter \x and on the arbitrary function /. 
In fact, it has the form of the following product 

j= (i+^'^iv*) (i+^'^ivA, (14) 



with qk given by (0). Notice that J takes values in SU(2); is smooth everywhere on R 2+1 
(mainly because, the two vectors qi and q2 are nowhere zero); it satisfies the boundary 
condition (H); and the equation of motion ([[]). 

To start with, and in order to illustrate the above family of soliton solutions, let us 
examine two simple cases, by giving specific values to the parameters \x and f(u). [The 
complex parameter /i determines the velocity of the "centre-of-mass" of the system.] 



Let us take fi = i (which corresponds to the "centre-of-mass" of the system being 
stationary) and f(u) = u>, thus u = z and <p — t, where z = x + iy; r 2 = zz. 
Therefore the row vectors (|T3D, become 

gi = (i+r 2 )(l, z) -2it{z,-l), 

fc = (M). Uj 
In this time-dependent solution, for t negative, a ring structure with reducing radius 
is obtained, which deforms to a single peak at t — and thereafter expands again to 
a ring. Figure 1 presents few pictures of the corresponding energy density at some 
representative values of time. Ring structures occur in the soliton scattering of many 
nonintegrable planar systems ||, H and are an approximation of two solitons. 

This picture can be confirmed by looking at the energy density of the solution, which 
is 

c _^ r 4 + 2r 2 + 4t 2 (2r 2 + l) + l 

[r* + 2r 2 + At 2 + l] 2 ' 1 } 

Notice that the energy density is time-reversible and rotationally symmetric (see 
below). For large (positive) t, the height of the ring (maximum of 8) is proportional 
to 1/t, while its radius is proportional to \fi. 

Accordingly, let us take /i = % and f(u) = uj 2 . Thus, the row vectors (|l^) are 

gi = (l+r 4 )(M 2 )-4^ 2 ,-l), 

q2 = (l,z 2 ). (17) 

Here, for negative t, a single peak occurs with an additional ring, which changes 
to a ring structure at t = and reverts back to the original form, for positive t 
(see Figure 2). However, these rings are not radiation since they travel with speed 
less than that of light. In fact, for large (positive) t, their velocity is approximately 
proportional to t~ 2 ^ 3 . [Note that we have set the velocity of the light, c, equal to 
the unity, so that in all our calculations we can use dimensionless quantities.] 



This leads to an energy density, which is 

r 10 + 18t 2 r 8 + 2r 6 + 4t 2 r 4 + r 2 + 2t 2 



£ = 64- 



[r 8 + 2r 4 + 16t 2 r 2 + l] 2 



Again, £ has the same symmetries as in ([16]). For large (positive) t, the height of 
the soliton peak is proportional to t 2 and its radius is proportional to 1/t; while the 
soliton ring spread out, becoming broader and broader, with height proportional to 
t -2 / 3 and radius proportional to t 1 ^ 3 . 



Finally, a general concluding remark should be made. Although ([I]) is not rotationally 



symmetric in the iq/-plane; when f(z) = z p the field J flT3|, 14) is invariant under the 
transformation z — > e^z, since 

/e*/ 2 \ /e-*/ 2 \ 
J J " 1 o e -^p/2 i I e*/ 2 i ' ^ ' 

This transformation does not affect the equation of motion ([j]) due to the chiral symmetry 
J — > /«Jr where k and r are constant SU(2) matrices. The main features of this time- 
dependent solution may be inferred as follow. If r is large, the field J is close to its 
asymptotic value J , as long as 2t/'/|/| 2 — > 0. But as 2t|/'|/|/| 2 rs 1, J departs from its 
asymptotic value Jo and a ring structure emerge with radius proportional to (2tp) 1 ^ p+1 \ 

III. Soliton-Soliton Scattering. 

We now move on to the more interesting question of scattering processes. In fact, we will 
use the method of section 2 to construct solutions of ([!]) representing scattering solitons. 
We will see that, in all head-on collisions of iV moving solitons the scattering angle is 
tt/N. Moreover, when the iV solitons are very close together, and in particular, when 
they are on top of each other the N lumps which represent them merge together to form 
a ring-like structure. Then, instead of moving towards the centre, they emerge from 
the ring in a direction that bisects the angle formed by the incoming ones. As we have 
already mentioned this nontrivial scattering is not usual in an integrable theory, but is 
exceptional. 

The scattering solutions arise if we take a solution of the simple-pole case (§) with n — 2, 
put Hi = \i + e, fi2 = H — £ and take the limit e — ► 0. The constraint 72(^2) — /i(^i) - ► 
as e — > has to be imposed, in order for the resulting solution ip to be smooth for all 



(u, v ,x). So let us write /i(a>i) = /(^i) + £h(u>i), ^2(^2) = /(<^2) — £^(^2), where / and 
/i are both rational functions of one variable (the examples of the previous section had 
h = 0). Once again J is given by flT3]), with the two- vectors q k given by 

ft = (l + |/| 2 )(l,/)+<?0i -/*)(/,-!), (20) 

?2=(1,/), 

where 1? = + /i(u;). So this solution belongs to a large family, since one may take / and 
h to be any rational meromorphic functions of to. Note that J is smooth on R 2+1 and 
satisfies its boundary condition, irrespective of the choice of / and h. 

It may seem strange that one can take the limit of a family of soliton solutions with 
trivial scattering, and obtain a new one with nontrivial scattering. Thus, it is interesting 
to study how the solitons are affected by varying e. To do so, let us take a solution of the 
simple-pole case (H) with n = 2, put fii = i + e, fi 2 = i — £, while taking f k = u; fc ; and 
study how the configuration of the two initial well separated solitons changes as e — > 
at a fixed time (t = —15). Figure 3 shows that as e — > the solitons disperse, shift and 
interact with each other. In other words, their internal degrees of freedom as well as the 
impact parameter change in this limit, making the process highly nontrivial. 

As an example, let us present two typical cases. 



• Let us take // = i, f(u) = uj and h(u) = u 3 ; thus 1? = t + z 3 . For r large, J is equal 
to its asymptotic value J , as long as djz 3 = 1 + t/z 3 w 1, but as z approaches 
any of the three cube roots of —t, then d — > 0, while J departs from J , and three 
localized solitons emerge. For t negative, the three solitons are approximately at the 
points: ((— t) 1 / 3 ^), (— (— t) 1 / 3 , ±y/3 (— t) 1 / 3 ) ; while for t positive, the solitons are 
at (-^,0), (f/ 3 ,±V3tV 3 ). 

More information can be deduced from the energy density, which is 

£= 16[2r 8 + 16r 6 + 19r 4 + 2r 2 (l + 8xy 2 t)+4t 2 (l + 2r 2 ) + l + 8x?/ 4 t 

-8xH- 16tx(x 2 -y 2 )]/[4r 6 + r 4 + 2r 2 + 4t 2 + 1 + 8tx(x 2 -3y 2 )} 2 . ^ ' 

The density £ is symmetric under the interchange t > —t, x 1— » —a; and y 1— > —y. 
For small (negative) t, the solitons form an intermediate state having the shape of 
a ring with three maxima on the direction of the incoming solitons which deforms 
to a circularly-symmetric ring at t — and then energy seems to flow around, until 
three other maxima are formed in the transverse direction, for small (positive) t. 



Figure 4 shows clearly the intermediate states with three maxima. The three new 
maxima then give rise to three new solitons emerging at 60° to the original direction 
of motion. During the intermediate phase solitons lose their identity. 

Finally something has to be said about their size. For large (positive) t, their height 
is proportional to t~ 4 / 3 , their radius is proportional to t 1 / 3 , while their speed is 
proportional to t~ 2//3 : therefore, they spread out and slow down. 

• Accordingly let us take fj, — i while choose f(cu) = u 2 and h{uj) = u 3 . Here J 
departs from its asymptotic value J , when z approaches the values ±y/—2t or zero 
(since d = z(2t + z 2 ) — > 0); and (again) three localized solitons emerge. In this case 
though, if t is negative, all three of them are on the x-axis at x ss ±y/—2t and at the 
origin; while if t is positive, they are on the y-axis at y ~ ±\/2t and at the origin. 
So the picture consists of three solitons: a static one at the origin, with the other 
two accelerating towards the origin, scattering at right angles and then decelerating 
as they separate. 

This can be observed from the energy density, which is 

£ = 32[r 12 + 2r 2 (r 8 + r 6 + 1) + 36t 2 r 8 + 4r 6 + 9r 4 + 8t 2 r 4 + 4t 2 + 12t(x w - y 10 ) 
+4t{x 2 - y 2 ){3 + 2x 2 y 2 + 6x 4 y A ) + 4t{x 6 - y 6 ){9x 2 y 2 - 2) - y 10 )]/[r 8 + 4r 6 
+2r 4 + 16tr 2 (t + x 2 - y 2 ) + l] 2 . 

(22) 

Here S is symmetric under the interchange t h- > —t, x # y\ therefore the collision 
is time symmetric, with the only effect the 90° scattering (no phase shift; no radi- 
ation). For large (positive) t, the height of the static soliton is proportional to t 2 
and its radius is proportional to 1/t; while the moving solitons expand with height 
proportional to t~ 2 / 3 and radius proportional to t 1 / 3 . 

In Figure 5 we present some pictures of the total energy densities of three solitons 
during a typical nontrivial evolution. 

In principle one should be able to visualize the emerging soliton structures when f{uj) = uj v 
and h{uj) = u g , i.e. are rational of degree p,g 6 N, respectively. In fact, for q > p the 
configuration consists of (p — 1) static solitons at the "centre-of-mass" of the system 
(if more than one, a ring structure is formed) accompanied byiV = g — p + 1 solitons 
accelerating towards the ones in the middle, scattering at an angle of n/N, and then 
decelerating as they separate. This follows from the fact that the field J departs from its 



asymptotic value Jo when $ = lu^- 1 )^ (u — H~ 2 v) + u N ) — > 0, which is true when either 
k>(p-i) = o or u;^ + p (u — ^~ 2 v) = 0; and this is approximately where the solitons are 
located. 

We conclude this section by investigating the corresponding case where ip(\) has a triple 
pole (and no others). Therefore, it is taken to have the form 

(23) 

As we have already mentioned, the reality condition (|5|) is satisfied if and only if ip 
factorizes into three simple factors of the following type 

V (a-m) Ikilr / V \m\\ J \ INI 2 / 

for some two- vectors ft. The requirement that the matrices A = (Lip)^ 1 and B = 
(Mip)^ -1 should be independent of A imposes differential equations on ft; which are 
three nonlinear equations, and it seems difficult to find their general solution. 

One way of proceeding is to take a solution for the simple-pole case (0) with n = 3, put 
fix = i + e, fi2 = i, yU3 = i — £ and take the limit e — > 0. In order to obtain a smooth 
solution -0 for all (u,v,x), it is necessary that — ^2(^2) — * 0, /i(k>i) — ^3(^3) — > 0, 

^2(^2) - ^3(^3) -> as e -> 0. So let us write = /(wi) + e/i^) + £ 2 #(o;i), 

72(^2) = fM, /sM = fM - eh(u a ) + £ 2 #(^ 3 ), where /, h and 5- are rational 
functions of one variable. On taking the limit, we obtain a ip of the form (|24"D , smooth on 
R 2+1 and such that the matrices A and B be independent of A. 

Consequently, J = 0(O) _1 is a smooth solution of (HI) of the form 

(25) 

with ft being in terms of f{z), h(z) and g(z) by 

ft = (1 + |/| 2 ) 2 (1, /) - 4z(6 + id){\ + |/| 2 )(/, -1) - 46 2 (f, -/ - 2z6) - 8id5(l, /), 
g 2 = (l + |/| 2 )(l,/)-2^(/,-l), 
93 = (1,/), 

(26) 

where b = tf'(z) + h(z) and d = t 2 f" (z) /2 + i(t - ?/)/' (2) /2 + th'(z) + g(z). Note that the 
two- vectors q2, ft here correspond to the ones given by (ED) for /i = i, respectively. 



Let us examine a sample example of this solution, since we may take /, h and g to be 
any rational meromorphic function of z. 

1 1 



• Let us take f(z) = 0, h(z) = z and g(z) = z 2 ; thus b = z and d = t + z 2 . This 
solution consists of two solitons coming in along the y-axis merging to form a peak at 
the origin and then two new solitons emerging along the x-axis. Figure 6 illustrates 
what happens near t — 0. 

The energy density of the system is, 

[32r 4 + 12r 2 - lQyr 2 + lQt 2 + lQty + 32t(x 2 - y 2 ) + l] 2 ' ^ ' 
which has a reflection symmetry around the x-axis. For large (positive) t, S is 
peaked at two points on the y-axis, namely y ~ Moreover, the height of the 

corresponding solitons is proportional to 1/t, and their radius is proportional to \/t; 
which means that the y-axis asymmetry vanishes at t —>■ oo. 

IV. Construction of Soliton-Antisoliton Solutions. 

In this section we construct a large family of solutions which as we will argue later, can 
be though of as representing soliton-antisoliton field configurations. Roughly speaking, 
solitons correspond to / being a function of the variable z, and antisolitons correspond 
to a function of z. 

One way to generate a soliton-antisoliton solution of (|I|), is to assume that has the 
form 

, n 1 ® m 1 n 2 ®m 2 . . 

*m = 1 + 71377 + 7A+i)"' (28) 

Here m fc for fc = 1, 2 are complex-valued two- vector functions of (t, z, z) (not depending 
on A). 

The idea is to find the n\, m\, ... such that the reality condition (||) holds, and such 
that the matrices A = (L?/>)-?/> _1 and = (Mip)ijj~ l are independent of A. One way of 
proceeding is to take the solution @ with n = 2, put /ii = z + e, //2 = — * — £ and take 
the limit e —>■ 0. In order for the resulting -0 to be smooth on R 2+1 it is necessary to 
take fi = fiyji), f 2 = — 1/ f(uj 2 ) — sh(uj 2 ), where / and h are rational functions of one 
variable. On taking the limit e — * 0, we then obtain a ^ as in (|28|) with m k = (m^m^) 
being holomorphic functions of z (or z), through the relations m 1 = (1, /), m 2 = (— /, 1), 
while 

2^(1 + |/| 2 ) _! , 2^; 
(1 + |/| 2 ) 2 +| W | 2 ^ (1 + |/| 2 ) 2 +|«;| 2 ' 



J 



;i+l/i 2 ) 2 +H 



, (31) 



2 _ 2 W _! 2^(1 + |/| 2 ) 2 

" " (1 + 1/| 2 ) 2 + |ti;| 2 "* (l + l/IT + M 2 ™ ' 1 J 

with 

iu = h/ 2 + 2tf. (30) 

So we generate a solution J = = 0) _1 , which depends on the two arbitrary rational 
functions / = f(z) and h = h(z). This solution has the form 

M 2 + 2i(/w + /w) - (1 + |/| 2 ) 2 - 2i (w - / 2 w) 
-2i (w - / 2 w) |uf - 2%{fw + /w) - (1 + |/| 2 ) 2 

with w given by (|30|) . In general, by taking f(z) = z p and h(z) = z q where p is a 
positive integer and q is a non- negative integer; the energy, obtained by integrating fl2|), 
is -E 1 = (2p + q)8ir. Roughly speaking, the solution looks like (2p + q) lumps at arbitrary 
positions in the xy-pl≠ which as we are going to see are a combination of solitons and 
antisolitons. 

A topological charge may be defined for the field J (|3lD by exploiting the connection of 
it with the 0(3) cr-model. The unmodified chiral model [i.e., (1) with V a = (0,0,0)] is 



equivalent to the 0(4) cr-model p0| through the relation 

J = /0o + *cr-0, (32) 

where cr are the usual Pauli matrices, and ((f) ,(j>) = (<f) , 0i, 2 , 3 ) is a four vector of 
real fields that are constrained to lie on S 3 , i.e. 4>q + <fi ' 4> = 1- The only static finite 
energy solutions of the 0(4) cr-model correspond to the embedding of the 0(3) cr-model 
21| . Therefore the only static solutions of (1) are the 0(3) embeddings that we shall 



describe. This is because for the one-soliton solution (static or Lorentz boosted in the 
y-axis) the term in (|1|) proportional to V a is zero, so the system behaves like the 0(4) 
model, for which the 0(3) embedding is totally geodesic. [However, for time- dependent 
configurations, the term proportional to V a is non-zero and will affect the evolution of 
the field, which will in general not lie in an 0(3) subspace of 0(4).] 

To proceed further, let us mention the topological aspects of the 0(3) and 0(4) a- models. 
In studying soliton-like solutions, we require that the field configuration has finite energy. 
This implies that the field must take the same value at all points of spatial infinity, so 
that space is compactified from R 2 to S 2 . At fixed time, the field is a map from S 2 into 



the target space. Now for the 0(3) model, the field is a map <f> '■ S 2 — > S 2 , and due to the 
homotopy relation 

vr 2 (S 2 ) = Z, (33) 

such maps are clasified by an integer winding number M which is a conserved topological 
charge. An expression for this charge is given by 

AT = (Syr)" 1 J ey 4> ■ (d i( f> A dj4>) d 2 x, (34) 

where i — 1, 2 with a;* = (x, y). 

Although, for the 0(4) model [the same argument is valid for ([I]) due to the topolog- 
ical aspects of the theory] the field at fixed time is a map (0o, <f>) '■ S 2 — > S 3 and the 
corresponding homotopy relation is 

vr 2 (S 3 ) = 0, (35) 

so there is no winding number. However, for soliton solutions that correspond to some 
initial embedding of 0(3) space into 0(4), there is a useful topological quantity, as we are 
going to see. 

Consider the 0(4) configuration which at some time corresponds to an 0(3) embedding, 
which we choose to be 0o — for definiteness. At this time the field is restricted to an 
S 2 equator of the possible S 3 target space. Suppose that the field never maps to the 
anti-podal points {^4i,^4 2 } = {</>o = 1> 0o — — 1} & t an y time, so the target space is 
Sq = S 3 — {Ai, *4. 2 }. Now Sq ~ S 2 x R, and thus we have the homotopy relation 

vr 2 (S 3 ) = tt 2 (S 2 X R) = tt 2 (S 2 ) © 7T 2 (R) = Z, (36) 

and therefore a topological winding number exists. An expression for this winding number 
is easy to give, since it is the winding number of the map after projection onto the chosen 
S 2 equator, i.e. 

W = (87T)- 1 j eij cj>' ■ {difi A djtf/) d 2 x, (37) 

where 4>' = 4>/\4>\- If the field does map to the anti-podal points {^li,^4 2 } at some time 
the winding number is ill defined at this time and if considered as a function of time M' 
will be integer valued but may suffer discontinuous jumps as the field moves through the 
anti-podal points. In the following examples, before comparing the solution J given by 
(j3l"l) with the 0(3) embedding it is convenient to perform the transformation J — » M J 



with M = (V2)- 1 
embedding. 



1 1 
-1 1 



so that the evolution of the field remains close to the 0(3) 



V. Soliton-Antisoliton Scattering. 

Usually in the nonintegrable models, there is an attractive force between solitons of op- 
posite topological charge. In fact, if the solitons and antisolitons are well separated, then 
they attract each other and eventually annihilate into a wave of pure radiation which 
spreads with the velocity of light |3], || . However, the interaction forces between solitons 
and antisolitons do depend on their configuration; in particular, they depend on the rel- 
ative orientation between them in the internal space. Therefore, the cross section for the 
soliton-antisoliton elastic scattering is non-zero. [In the real world, the proton-antiproton 
elastic scattering is seen in a reasonable fraction of the cases.] This is the first example 
for which there has been constructed an explicit (since the system is integrable) solution 
of elastic soliton-antisoliton scattering in either integrable or nonintegrable model. As a 
result, it provides a major link between soliton dynamics in integrable and nonintegrable 
systems. 

The evolution is initially similar to the numerical results obtained through the connection 
of the integrable chiral model ([[]) with the 0(3) cr-model [2U|. In particular, a soliton 



and an antisoliton are moving along the x-axis towards each other at an accelerating 
rate until they merge at the origin and form a peak. Note that a peak is formed rather 
than a ring since the energy is mainly kinetic when a soliton and an antisoliton merge. 
However, rather than the peak dissipating in a wave of radiation it now reforms into two 
new structures which undergo 90° scattering. In general, in all head-on collisions of iV 
moving soliton and antisoliton objects, the scattering angle is ir/N degrees relative to the 
initial direction of motion. 

Next we looked at two cases corresponding to the mixtures of solitons and antisolitons. 
[The configurations given by (|3l|) when h(z) = are equivalent to the ones obtained from 
when f(z) = z p .} 



• First, let us take f(z) = z and h(z) = 1. Roughly speaking, if r is large, J is close 
to its asymptotic value J , as long as w/z 2 = 1 + 2t/z 2 ps 1; but as z approaches 

1 rr 



±y/—2t then w — > 0, and J departs from its asymptotic value: this is where the 
two structures are located. More precisely, for negative t, the two objects are on the 
x-axis, approximately at x ~ ±\/—2t; while for positive t, they are on the y-axis, 
approximately at y ~ ±\/2£. Figure 7 illustrates what happens near t = 0. 

The picture is consistent with the properties of the energy density of the solution, 
which is 

2r 4 + 4r 2 + 4t 2 (l + 2r 2 ) - 4t(x 2 - y 2 ) + 1 

[2r 4 + 2r 2 + At{x 2 - y 2 ) + At 2 + l] 2 ' { ' 

Note the symmetry of £ under the interchange t i— > — t, x ^ y; the time symmetry 

of the density confirms the lack of radiation. The corresponding localized structures 

are not however of constant size: for large (positive) t, their height is proportional 

to 1/t, while their radius is proportional to \rt. 

The projected topological charge M' is zero throughout the scattering process; while 
the projected topological density q', i.e. 

W = J q dx dy, (39) 

has an almost identical distribution (up to a scale) to that of the energy density (see 
Figure 8(a)). Therefore, the configuration represents a soliton and an antisoliton 
which are clearly visible as distinct structures having respectively +1 and -1 units 
of topological charge concentrated in a singe lump. 



Equation ([!]) is not Lorentz invariant and indeed is not even radially symmetric due to the 
presence of the vector V a which picks out a particular direction in space, and therefore 
one may expect to find different scattering behaviour for more general solutions; e.g., 
when the soliton and the antisoliton are moving along the x-axis rather that the y-axis. 
However, this is not true since ([!]) is a reduction of the self-dual Yang-Mills equation 
in R 2+2 which does have an SO(l,2) symmetry. Therefore, the SO(2) symmetry of the 
Yang-Mills system means that any given solution J, can in principle be converted to 
gauge fields by performing a coordinate rotation (together with a gauge transformation) 
and then recover the corresponding J' which will describe the same solution as J but with 
a rotated coordinate system. Indeed, this is what happens by taking 

f( z ) = e ^z, h{z) = 1, (40) 

where is an angle in the xj/-plane. This picture presents a rotated version through any 
angle <fi in the x?/-plane of the original one (i.e., Figure 7). 



• Finally, let us take f(z) = z and h(z) = z. The corresponding configuration consists 
of one soliton and two antisolitons (see Figure 8(b)). 

It is interesting to look at the time dependence of various energies in each process. 
The total energy, of course, is constant and it is the spatial integral of the following 
energy density 

£ = 8[r 8 + 8r 6 + llr 4 + 4r 2 - 8xH + 16ty 2 (x 3 + t) + 8t 2 + A8xy 2 t + 2 

-16x 2 t(x-t) + 24xty*}/[r 6 + r 4 + 2r 2 + 4t 2 + 4tx 3 - 12xy 2 t + l] 2 . ^ ^ 

Obviously, the energy density £ is symmetric under the interchange t ^ —t, x i— > —x 

and y i— > — y, only. Again all three structures come together forming a bell-like 

structure and then emerge at an angle of 60° with respect to the original direction. 

However, by looking at the maximum of £ we observe that, for large (positive) t, 

the height of the localized structures is proportional to t~ 4 / 3 , while their radius is 

proportional to t 1//3 ; thus they spread out as they move apart. 

Figure 9 shows the results of a head-on collision of the one-soliton two-antisoliton 
system. 

Let us conclude with the observation that, by taking f(z) = z p and h(z) = z q , J departs 
from its asymptotic value Jq when w = z p ~ 1 (2tp + z ) — > with N = p + q + 1, which is 
true when either z^ 1 ^ = or 2tp + z N = 0: this is approximately where the lumps are 
located. Therefore, J represents a family of soliton-antisoliton solution which consists of 
(p — 1) static soliton-like objects at the origin, with iV others accelerating towards them, 
scattering at an angle of tc/N, and then decelerating as they separate. 

VI. Conclusion. 

The infinite number of conservation laws associated with a given integrable system place 
severe constraints upon possible soliton dynamics. The construction of exact analytic 
multisoliton solutions with trivial scattering properties is a result of such integrability 
properties. In this paper new soliton and soliton-antisoliton solutions have been obtained 
for the planar integrable chiral model (|l]). These structures travel with non-constant 
velocity; their size is non-constant; and they interact non-trivially. Such results might be 
useful for connecting integrable and nonintegrable systems which possess soliton solutions. 
In addition, they indicate the likely occurrence of new phenomena in higher dimensional 
soliton theory that are not present in (1+1) dimensions. 



It seems likely that there are many more interesting solutions still to be found; an open 
question being what is the general form of the function ip when it has a higher-order pole 
in A. One could, for example, investigate the case n = 3 for with a single and a 
double pole; and determine the scattering properties of the emerging structures, in terms 
of their initial velocity and of the values of the impact parameter. Finally, it would be of 
great interest to deduce the general form of the function for the soliton-antisoliton 
case ([28]) with the only constraint to satisfy the reality condition (|5|) and the requirement 



that the matrices A = (Lip)^ 1 and B = (Mifj)^ -1 be independent of A. 
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Figure Captions. 

Figure 1: The energy density S (16) at increasing times. 
Figure 2: The energy density £ (18) at various times. 

Figure 3: Energy density at various values of e for a system of two solitons (t = —15). 

Figure 4: Energy density at increasing times for a system of three solitons with 60° angle 
scattering. 

Figure 5: Energy density at various times for the scattering of three solitons, with one 
being static at the origin. 

Figure 6: Energy density at increasing times when tp(X) has a triple pole (and no others). 

Figure 7: Energy density at increasing times showing a 90° scattering between a soliton 
and an antisoliton. 

Figure 8: Topological charge density at increasing times for (a) soliton-antisoliton scat- 
tering, (b) one-soliton two-antisoliton scattering. 

Figure 9: Energy density of a system consisting of a soliton and two antisolitons at 
various times. 
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